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Abstract 

A new class of LRS Bianchi type VIo cosmological models with free 
gravitational fields and a variable cosmological term is investigated in 
presence of perfect fluid as well as bulk viscous fluid. To get the deter- 
ministic solution we have imposed the two different conditions over the 
free gravitational fields. In first case we consider the free gravitational 
field as magnetic type whereas in second case 'gravitational wrench' of 
unit 'pitch" is supposed to be present in free gravitational field. The vis- 
cosity coefficient of bulk viscous fluid is assumed to be a power function 
of mass density. The cosmological constant A is found to be a decreasing 
function of time and positive which is corroborated by results from recent 
supernovae la observations. The physical and geometric aspects of the 
models are discussed. 

PACS: 98.80.Cq, 98.80.-k 
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1 Introduction and Motivations 

The problem of the cosmological constant is salient yet unsettled in cosmol- 
ogy. The smallness of the effective cosmological constant recently observed 
(Ao < 10~ 56 cm~ 2 ) poses the most difficult problems involving cosmology and 
elementary particle physics theory. To explain the striking cancellation between 
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the "bare" cosmological constant and the ordinary vacuum energy contribu- 
tions of the quantum fields, many mechanisms have been proposed pQ. The 
"cosmological constant problem" can be expressed as the discrepancy between 
the negligible value of A for the present universe as seen by the successes of 
Newton's theory of gravitation [2 whereas the values 10 50 larger is expected by 
the Glashow-Salam- Weinberg model [3J and by grand unified theory (GUT) it 
should be 10 107 larger [4]. The cosmological term A is then small at the present 
epoch simply because the universe is too old. The problem in this approach is 
to determine the right dependence of A upon R or t. 

Models with a relic cosmological constant A have received ample attention 
among researchers recently for various reasons (see Refs. [5] - [10] and references 
therein). Some of the recent discussions on the cosmological constant "prob- 
lem" and on cosmology with a time- varying cosmological constant by Ratra and 
Peebles [11], Dolgov et al. [12 HE], Dolgov [14], and Sahni and Starobinsky [15] 
point out that in the absence of any interaction with matter or radiation, the 
cosmological constant remains a "constant" , however, in the presence of interac- 
tions with matter or radiation, a solution of Einstein equations and the assumed 
equation of covariant conservation of stress-energy with a time-varying A can 
be found. For these solutions, conservation of energy requires decrease in the 
energy density of the vacuum component to be compensated by a correspond- 
ing increase in the energy density of matter or radiation. Earlier researches 
on this topic, are contained in Zeldovich [16], Weinberg [1] and Carroll, Press 
and Turner [17] • Recent observations by Perlmutter et al. [18] and Riess et 
al. [19] strongly favour a significant and positive value of A. Their findings 
arising from the study of more than 50 type la supernovae with redshifts in the 
range 0.10 < z < 0.83 suggest Friedmann models with negative pressure matter 
such as a cosmological constant (A), domain walls or cosmic strings (Vilenkin 
[20] . Garnavich et al. [21]). Recently, Carmeli and Kuzmenko [22] . Behar and 
Carmeli [53] have shown that the cosmological relativistic theory predicts the 
value for cosmological constant A = 1.934 x 10~ 35 s~ 2 . This value of "A" is 
in excellent agreement with the measurements recently obtained by the High-Z 
Supernova Team and Supernova Cosmological Project (Garnavich et al. [2"T] , 
Perlmutter et al. [18 , Riess et al. [19], Schmidt et al. [24]) The main conclusion 
of these observations is that the expansion of the universe is accelerating. 

Several ansiitz have been proposed in which the A term decays with time 
(see Refs. Gasperini [25] , Berman [3DJ, Berman et al. [17], Ozer and Taha [7J, 
Freese et al. [8], Peebles and Ratra [28], Chen and Wu [29], Abdussattar and 
Viswakarma [3D], Gariel and Le Denmat [3T], Pradhan [32, Pradhan et al. [33]). 
Of the special interest is the ansatz A cx S~ 2 (where S is the scale factor of the 
Robertson- Walker metric) by Chen and Wu [25], which has been considered or 
modified by several authors (Abdel-Rahaman [33], Carvalho et al. [5J, Silveira 
and Waga [10], Vishwakarma [35]). 

In classical electromagnetic theory, the electromagnetic field has two inde- 
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pendent invariants FijF^ and *FijF l K The classification of the field is char- 
acterized by the property of the scalar k 2 = {F^F'^) 2 + (* F^F^) 2 ■ When k 
= the field is said to be null and for any observer \E\ = \H\ and E.H = 
where E and H are the electric and magnetic vectors respectively. When k 7^ 0, 
the field is non-null and there exists an observer for which mE = nH, m and 
n being scalars. It is easy to see that if = 0, FijF^ ^ then either 

E = or H = 0. These we call the magnetic and electric fields, respectively. 
If F i:j F l 3 = QfFijF** ^ then E = ±H. This we call the 'electromagnetic 
wrench' with unit 'pitch'. In this case Maxwell's equations lead to an empty 
electromagnetic field with constant electric and magnetic intensities. In the case 
of gravitational field, the number of independent scalar invariants of the second 
order is fourteen. The independent scalar invariants formed from the conformal 
curvature tensor are four in number. In the case of Petrov type D space-times, 
the number of independent scalar invariants are only two, viz. ChijkC hl i k and 
*ChijkC hl; > k . Analogous to the electromagnetic case, the electric and magnetic 
parts of free gravitational field for an observer with velocity if are defined as 
E a f3 — C a jpiV l v^ and H a p = *C a jpiV l v^ ( [36J ) . It is clear from the canonical 
form of the conformal curvature for a general type D space-time that there exists 
an observer for which E a p = {—)H a p, where m,n being integers and m 0. 
The field is said to be purely magnetic type for n = 0, m ^ 0. In this case we 
have E a p = and H a p ^ 0. The physical significance for the gravitational field 
of being magnetic type is that the matter particles do not experience the tidal 
force. When m ^ and also n ^ 0, we call that there is a 'gravitational wrench' 
of unit 'pitch' | — | in the free gravitational field [37] . If 'pitch' is unity then we 
have E a p = ±H a/3 . 

The space-time having a symmetry property is invariant under a continuous 
group of transformations. The transformation equations for such a group of 
order r is given by 

= P{x\...,x r ,a\...,a r ) (1) 
which satisfy the differential equations 

— =C {0) {X)Ai{a), (a,0 = l,...,r) (2) 

where a 1 , a r are r essential parameters. The vectors are the Killing vectors 
for the group G r of isometry satisfying the Killing's equation 

£(a)»;j + = ( 3 ) 

A subspace of space-time is said to be the surface of transitivity of the group 
if any point of this space can be transformed to another point of it by the 
action of this group. A space-time is said to be spatially homogeneous if it ad- 
mit a group G r of isometry which is transitive on three dimensional space-like 
hyper-surfaces. The group G3 of isometry was first considered by [38] who ob- 
tained nine different types of isometry group known as the Bianchi types. The 
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space-time which admits G4 group of isometry is known as locally rotationally 
symmetric (LRS) which always has a G3 as its subgroup belonging to one of 
the Bianchi type provided this G3 is simply transitive on the three dimensional 
hyper-surface t = constant. 

Considerable work has been done in obtaining various Bianchi type cosmo- 
logical models and their inhomogeneous generalization. Barrow [39j pointed 
out that Bianchi VIo models of the universe give a better explanation of some 
of the cosmological problems like primordial helium abundance and they also 
isotropize in a special sense. Looking to the importance of Bianchi type VIo 
universes, many authors [10] — [H] have studied it in different context. Recently 
Bali et al. [45] have obtained some LRS Bianchi type VIo cosmological models 
imposing two types of conditions over the free gravitational fields. 

In this paper we have revisited and generalized the solutions obtained by 
Bali et al. [JS]. We have considered an LRS Bianchi type VIo space-time and 
obtained models with free gravitational field of purely 'magnetic type' and also 
in the presence of 'gravitational wrench' of unit 'pitch' in the free gravitational 
field. It is found that the 'magnetic' part of the free gravitational field induces 
shear in the fluid flow, which is zero in the case of a 'electric' type free grav- 
itational field representing an unrealistic distribution in this case. This paper 
is organized as follows. The introduction and motivation are laid down in Sec. 
1. The metric and the field equations are given in Sec. 2. In Sec. 3, solutions 
representing LRS Bianchi type VIo cosmological models with perfect fluid and 
bulk viscous fluid are obtained imposing the condition when the free gravita- 
tional field is purely magnetic type (m 7^ 0, n = 0. In Sec. 4, we obtain the 
solution in presence of perfect fluid imposing the condition when there is a 'grav- 
itational wrench' of unit 'pitch' in the free gravitational field i.e. E a p — ±H a p. 
Discussion and concluding remarks are given in the last Sec 5, . 

2 The Metric and Field Equations 

We consider an LRS Bianchi type VIq universe for which 



where 9 1 = A(t)dx, 6 2 = B(t) exp (x)dy, 9 3 = B(t) exp (-x)dz, 6> 4 = dt. 

The energy-momentum tensor for a perfect fluid distribution with comoving 

flow vector v l is given by 



where v l — 8\, p and p being respectively, energy density and thermodynamic 
pressure of the fluid. Here we obtain 




(4) 



Tj = (p + p)viV j + pSj 



i ' 



(5) 



T 1 = T$ = Tf = p, T 4 4 



= -P- 



(6) 



4 



The Einstein's field equations (in gravitational units c = 1, G = 1) read as 
Ri-lR6i+A5i = -8nT?, (7) 
for the line element (H|) has been set up as 

A = -8vrp, (8) 

+ A = -8np, (9) 
+ A = 8np. (10) 
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Here, and also in the following expressions a dot indicates ordinary differentia- 
tion with respect to t. 

The average scale factor S for LRS Bianchi type VIo model is defined by 

S=(AB 2 )i. (11) 

A volume scale factor is given by 

V = S 3 = (AB 2 ). (12) 
The generalized mean Hubble parameter H is given by 

H = ^(H 1 +H 2 + H 3 ), (13) 
where i?! = ^ H 2 = H 3 = § 

The expansion scalar 9 and shear scalar er are obtained as 



and 



(14) 



respectively. The average anisotropy parameter is given by 



1 — 1 N 



where AHt = Hi - H(i = 1, 2, 3). 

The deceleration parameter (q) is defined as 



s 
s 

s 2 ' 



(17) 
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The non-vanishing physical components of Cijki for the line (|18p are given by 
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Equations ([8|)- p0|l are three relations in five unknowns A, B, p, p and A. 
For complete solutions of equations (|5|)- (fTU)) . we need two extra conditions. The 
research on exact solutions is based on some physically reasonable restrictions 
used to simplify the Einstein equations. Although the distribution of matter 
at each point determines the nature of expansion in the model, the latter is 
also affected by the free gravitational field through it's effect on the expansion, 
vorticity and shear in the fluid flow. A prescription of such a field may therefore 
be made on an a priori basis. The cosmological models of Friedman Robertson 
Walker, as well as the universe of Einstein - de Sitter, have vanishing free grav- 
itational fields. In the following two cases we impose different conditions over 
the free gravitational field to find the deterministic solutions. 



3 First Case: Free Gravitational Field is Purely 
Magnetic Type 

3.1 Solution in Presence of Perfect Fluid 

In this section we have extended the solution obtained by Bali et al. [45] by 
revisiting their solution. When free gravitational field is purely magnetic type 
(m ^ 0, n = 0), we have H a p ^ and E a p = 0. From (fT8|) . we obtain 



.4 
A 



B 
B 



AB 
AB 



B 2 



A 2 



Equation 



together with © reduce to 
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From Eqs. (|2U)l and (|2"Tj) , we obtain two independent equations 



Let % = U. Then Eqs. and {23]) take the form 

UU& + UU^A - 2Uf = 0, 

and 



(24) 
(25) 



(26) 
(27) 



C/| - UU ( A + 2U 2 = 0, 
respectively, where £ is defined by 

d _ 1 d 
dt ~ Irf^' 

Eqs. ([M]) and (J^SJ) lead to 

UU# - C/| + 2C/ 2 = 0. 

On substituting [/ = exp(/i), the above second order differential equation re- 
duces to the form 

M« + 2 = 0, (28) 

which gives 

LL = k 1 t-e + k 2 , (29) 
where k\ and ki are arbitrary constants. Hence, we obtain 

U = k 3 cxp[-T(T~k 1 )}, (30) 

where T stands for £. Therefore, from Eqs. ([25]). (|25|) and ([30]) . we obtain 

fc 4 exp [-T(T-fci)] 



B = 



(fci - 2T) 
(fci - 2T) ' 



where T is given by 



dT _ 1 
dt ki 



(fci -2T)exp [T(r-fci)], 



(31) 
(32) 

(33) 



^3, fc 4 and fcs(= ||) being arbitrary constants. Therefore the geometry of the 
universe Q reduces to the form 



ds 



2 _ /fcfcxp [-2T(T- fci)] 
(fci - 2T)2 



■ rfT2 + da; 2 + exp[2T(r-fc 1 



{exp (2x)dy 2 + exp (-2x)d^ 2 } 



(34) 



The expressions for pressure p and density p for the model (|34|) are given by 
3 



87175 : 



[4 - (fci - 2T) 2 ] exp [2T(T - fci)] - A(T), 



(35) 
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8np = _ [4 + (fcx - 2T) 2 ] exp [2T(T - h)] + A(T), 



(36) 



For the specification of A(T), we assume that the fluid obeys an equation of 
state of the form 

P = IP, (37) 

where 7(o<7<l)isa constant. 
Using ([57)1 in ([33)) and (|3"6")l , we obtain 



8tt(1 + i)p = £(fc x - 2T) 2 cxp [2T(T - fc x ) 

^4 



(38) 



Eliminating p between Eqs. (1331) and ([33]) . we obtain 



(| - )A = — 12(1 + 7) exp [2T(T - fa)] + 3(1 a 7) (h - 2T) 2 exp [2T(T - h) 

Ka 



(39) 




Figure 1: The plot of energy density p(t) Vs. time 



From Eq. (|38[> . we observe that p(t) is a decreasing function of time and 
p > for all times. Figure 1 shows this behaviour of energy density. From Eq. 
(|39|) . we note the cosmological term A is a decreasing function of time and it 
approaches a small positive value with increase in time. From Figure 2 (pandA 
are in geometrical units in entire paper), we note the same character of A. This 
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Figure 2: The plot of cosmological term A(t) Vs. time 



is to be taken as a representative case of physical viability of the model. 

The behaviour of the universe in this model will be determined by the cos- 
mological term A, this term has the same effect as a uniform mass density 
p e ff = — A/4-7T which is constant in space and time. A positive value of A corre- 
sponds to a negative effective mass density (repulsion). Hence, we expect that 
in the universe with a positive value of A the expansion will tend to accelerate 
whereas in the universe with negative value of A the expansion will slow down, 
stop and reverse. In a universe with both matter and vacuum energy, there is a 
competition between the tendency of A to cause acceleration and the tendency 
of matter to cause deceleration with the ultimate fate of the universe depending 
on the precise amounts of each component. This continues to be true in the 
presence of spatial curvature, and with a nonzero cosmological constant it is no 
longer true that the negatively curved ("open") universes expand indefinitely 
while positively curved ( "closed" ) universes will necessarily recollapse - each of 
the four combinations of negative or positive curvature and eternal expansion 
or eventual recollapse become possible for appropriate values of the parame- 
ters. There may even be a delicate balance, in which the competition between 
matter and vacuum energy is needed drawn and the universe is static (non ex- 
panding). The search for such a solution was Einstein's original motivation for 
introducing the cosmological constant. Recent cosmological observations (Gar- 
navich et al. [21], Perlmutter et al. [18], Riess et al. [19], Schmidt et al. [24] ) 
suggest the existence of a positive cosmological constant A with the magnitude 
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A(Gh/c 3 ) « 10~ 123 . These observations on magnitude and red-shift of type la 
supernova suggest that our universe may be an accelerating one with induced 
cosmological density through the cosmological A-term. Thus the nature of A in 
our derived model of the universe is consistent with recent observations (Gar- 
navich et al. [21] . Perlmutter et al. [18] , Riess et al. [TI5], Schmidt et al. [21]). 



Some Physical and Geometric Features of the Model 

The expressions for kinematical parameters i. e. the scalar of expansion 9, 
shear scalar er, average scale factor S, proper volume V 3 and average anisotropy 
parameter A p for the model (|36|) are given by 

9=±-[6 + (k 1 - 2T) 2 ] exp [T(T - fei)], (40) 
' (k 1 -2T) 2 exp[T(T-k 1 )}, (41) 



5= (fc 1 -2T) 4 ex P ^™- ^ 



3 

, fc|expJ-22XT-fci)j 
V ~ V 5 ~ fc|(fc 1 -2T)4 ' 



A P = \- (44) 
The directional Hubble's parameters H x , H y and H z are given by 

H x = l-[(k! - 2T f + 2] exp [T(T - h], (45) 

H y =H Z = ^- exp [T(T - kx] , (46) 
where the mean Hubble's parameter is given by 

H =-^ r {6+(fc 1 -2T) 2 }exp[T(T-fc 1 )] (47) 

The model (|36p starts expansion with a big-bang singularity from T = — 00 and 
it goes on expanding till T = respectively correspond to the cosmic time 
t = and t = 00. The is found to be realistic everywhere in this time interval 
for A > — i§ exp (~~2^ • The model behaves like a steady-state de-Sitter type 
universe at late times where the physical and kinematic parameters p, p, 9 tend 
to a finite value, however shear vanishes there. The model has a point type 
singularity at time T = k\. The singular behaviour may be close to cosmic 
origin or outside the evolution. The average anisotropy parameter A p remains 
uniform and isotropic through out the evolution of the universe. This would 
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Figure 3: The plot of average scale factor S Vs. time 

depend on physical properties of matter and radiation. This may need detailed 
study to make better quantifiable view. From figure 3, it can be seen that in 
the early stages of the universe, i. e., t = 0, the scale factor of the universe had 
been approximately constant and had increased very slowly. At specific time 
the universe had exploded suddenly and expanded to large scale. This is good 
matching with big bang scenario. This is indicated in first part (top) of figure 
3. Later singular behaviour depends on (ki,T). 

3.2 Solutions For Bulk Viscous Fluid 

Astronomical observations of large-scale distribution of galaxies of our universe 
show that the distribution of matter can be satisfactorily described by a perfect 
fluid. But large entropy per baryon and the remarkable degree of isotropy of 
the cosmic microwave background radiation, suggest that we should analyze 
dissipative effects in cosmology. Further, there are several processes which are 
expected to give rise to viscous effect. These are the decoupling of neutrinos 
during the radiation era and the recombination era [46] , decay of massive super 
string modes into massless modes [47] , gravitational string production [48] [49] 
and particle creation effect in grand unification era [50] . It is known that the 
introduction of bulk viscosity can avoid the big bang singularity. Thus, we 
should consider the presence of a material distribution other than a perfect fluid 
to have realistic cosmological models (see Gr0n |51) for a review on cosmological 
models with bulk viscosity). A uniform cosmological model filled with fluid 
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which possesses pressure and second (bulk) viscosity was developed by Murphy 
[52] , The solutions that he found exhibit an interesting feature that the big 
bang type singularity appears in the infinite past. 

In presence of bulk viscous fluid distribution, we replace isotropic pressure p by 
effective pressure p in Eq. (|35p where 

P=p-£t£> (48) 
where £ is the coefficient of bulk viscosity. 

The expression for effective pressure p for the model Eq. (|34[) is given by 

8np = 8n(p - £<) = — J [4 - (k x - 2T) 2 ] exp [2T(T - fci)] - A(T). (49) 

Thus, for given £(t) we can solve for the cosmological parameters. In most of the 
investigation involving bulk viscosity is assumed to be a simple power function 
of the energy density (Pavon |53j . Maartens [54j . Zimdahl [55], Santos [56] ) 

f(t) = ZoP n , (50) 

where £o and are constants. For small density, n may even be equal to unity 
as used in Murphy's work [55] for simplicity. If n = 1, Eq. ([M]) may correspond 
to a radiative fluid (Weinberg .2 ). Near the big bang, < n < i is a more ap- 
propriate assumption (Belinskii and Khalatnikov |57p to obtain realistic models. 

For simplicity sake and for realistic models of physical importance, we con- 
sider the following two cases {n = 0, 1): 

3.2.1 Model I: When n = 

When n = 0, Eq. ([50]) reduces to £ = £o = constant. With the use of Eqs. ([3"6"]) , 
([37]) and (gOJ, Eq. gl]) reduces to 

8tt(1+ 7 )p= -^(fci -2T) 2 exp[2T(T-/ci)] + 

^{6 + (/e x -2T) 2 } exp [2T(T -**)]. (51) 
K4 

Eliminating p(t) between Eqs. ([3"6"]l and (|5 1 [) . we obtain 



(| , _ ) A : _ 12 (! + 7) exp [2T(r _ fci)] + 3(1 7) (fci _ 2T)2 exp [2y(T _ fci) 



^{6+(k 1 -2T) 2 }exp[T(T-k 1 )}. (52) 
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Figure 4: The plot of energy density p(t) Vs. time 




Figure 5: The plot of cosmological term A(t) Vs. time 
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Figure 6: The plot of energy density p(t) Vs. time 




Figure 7: The plot of cosmological term A(t) Vs. time 
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3.2.2 Model II: When n = 1 

When n = 1, Eq. flSOJ) reduces to £ = £ p. With the use of Eqs. (36} , J37} and 
|gDD , Eq. 105]) reduces to 



3(fca - 2T) 2 exp [2T(T - k x ) 



r 47rfc 4 [fc*(l + 7) - £o{6 + (fci - 2T) 2 } exp [T(T - fci)]] 
Eliminating pit) between Eqs. ([3"6]) and (|53[) . we obtain 

12 ro^v^ 3(fci - 2T) 2 exp[2T(T- fci)] 



(53) 



A =-T2 ex P[ 2T ( T - fc i)] 



fe 4 fc| 



fc 4 ( 7 + 1) - £o{6 + (fci - 2T) 2 } exp [T(T - h)} 



fc 4 (7 - 1) - £o{6 + (h - 2T) 2 } exp [T(T - fc x )]j (54) 
where fci, fc 4 , £oj and 7(0 < 7 < 1) are constants. 

From Eqs. (fSTj) and ((53^) . we observe that p(<) in both models are a decreas- 
ing function of time and p > for all times. Figures 4 and 6 show this behaviour 
of energy density. From Eqs. ([52]) and (f54|) . we see that the cosmological terms 
A in both models are a decreasing function of time and they approach a small 
positive value at late time. From Figures 5 and 7, we note this behaviour of 
cosmological term A in both models I and II. Thus, our models are consistent 
with the results of recent observations (Garnavich et al. [3T], Perlmutter et al. 
[15] . Riess et al. [19], Schmidt et al. [24]). 

The effect of bulk viscosity is to produce a change in perfect fluid and there- 
fore exhibits essential influence on the character of the solution. A comparative 
inspection of Figures 1, 4, 6 show apparent evolution of time due to perfect 
fluid and bulk viscous fluid. It is apparent that the vacuum energy density (p) 
decays much fast in later case. Also shows effect of uniform viscosity model 
and linear viscosity model. Even in these case decay of vacuum energy density 
is much faster than uniform. So the coupling parameter £ would be related 
with physical structure of the matter and provides mechanism to incorporate 
relevant property. In order to say more specific, detailed study would be needed 
which would be reported in future. Similar behaviour is observed for the cosmo- 
logical constant A. We also observe here that Murphy's [55] conclusion about 
the absence of a big bang type singularity in the infinite past in models with 
bulk viscous fluid in general, is not true. The results obtained by Myung and 
Cho [S5] also show that, it is not generally valid since for some cases big bang 
singularity occurs in finite past. For both models, it is observed that the effect 
of viscosity prevents the shear and the free gravitational field from withering 
away. 
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4 Solution of Field Equations for Second Case 

In this case there is a 'gravitational wrench' of unit 'pitch' in the free gravita- 
tional field i.e. E a p = ±H a p. Therefore we have 

E a/3 = KH af3 , k 2 = 1. (55) 

In this case, Bali et al. [IS] have investigated the solution given by 



ds 2 = 



(2t 2 + 3kt + 2)2 



3k 

exp I — = tan 
1 V7 



_! ( 4t + 3k) 
V7 



dT 2 



(2t 2 +3kt + 2)62 



/ .^tan' 1 (4T+ ^ K) ){e 2x dy 2 +e- 2x dz 2 } 

V7 V7 1 



-dx 2 + (2t 2 +3kt+2)5 exp 



(56) 



The expressions for pressure p and energy density p for the model (156|) are 
obtained as 

(12kt 3 + 29t 2 + 72kt-48) / 3k _ 1 (4t + 3k)\ ... 

8np = - ; exp = tan 1 - A(t) (57) 

4(2t 2 + 3kt + 2)^ V V7 y/7 J yl y 1 

(12kt 3 + 39r 2 + 72kt + 48) / 3k _ x (4t + 3k) \ A . . . . 

87rp = j exp = tan 1 = — - + AM (58) 

4(2t 2 + 3kt + 2)^ V V7 V7 J 

For the specification of A(r), we assume that the fluid obeys an equation of 
state of the form (137]). Using Eqs. ([37]) in ([57]) and 1(551), we obtain 



87r(l + 7)p: 



6kt 3 + 17r 2 + 36kt 



3k 

exp < = tan" 

1 V7 



4r + 3k 
~7T 



(2t 2 + 3kt + 2)2 
Eliminating p between ((57)) and ((59"|) . we obtain 

'12(7 ~ 1)«(t 2 + 6)t + T 2 (39 7 - 29) + 48( 7 + 1) 



(59) 



A = - 



4( 7 + 1)(2T 2 + 3kT + 2)5 



3k /4T + 3k 

exp < = tan = — 

1 V7 V V7 



(60) 



From preliminary study, we find that both energy density and cosmological 
constant are negative. Hence it will not be studied. It also shows singular 
behaviour in energy density at later stage of the evolution. Hence, the model is 
unphysical for further study. 



5 Discussion and Concluding Remarks 

In this paper, we have studied properties of the free gravitational field and their 
invariant characterizations and obtained LRS Bianchi type VIo cosmological 
models imposing different conditions on the free gravitational field. In first case, 
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where fee gravitational field is purely magnetic type, we observe that the energy 
density/? and cosmological constant A are well behaved. Also the effect of bulk 
viscous fluid distribution in the universe is compared with perfect fluid model. 
We observe that due to presence of bulk viscous fluid, the rate of decrease in 
energy density is faster compared to perfect fluid model. The linear relation 
of the coefficient of bulk viscosity with mass density provides further enhance 
decrease rate (see figures 1, 4, 6). Since we had considered extension of the 
various models, detailed physical parameter study would be reported in future. 
Important incorporation is time-dependence of cosmological constant A. The 
similar enhance decrease rate is also observed for A (see figures 2, 5, 7). The 
scale factor dependence is also reported. The nature of A in our derived models 
of the universe is found to be consistent with recent observations (Garnavich et 
al. [21], Perlmutter et al. [IS], Riess et al. QU, Schmidt et al. [24]). 
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